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ABSTRACT 

We present stndies of C/2015 D1 (SOHO), the hrst snnskirting comet ever 
seen from gronnd stations over the past half centnry. The Solar and Heliospheric 
Observatory (SOHO) witnessed its peculiar light curve with a huge dip followed 
by a flareup around perihelion: the dip was likely caused by sublimation of 
olivines, directly evidenced by a coincident temporary disappearance of the tail. 

The flareup likely reflects a disintegration event, which we suggest was triggered 
by intense thermal stress established within the nucleus interior. Photometric 
data reveal an increasingly dusty coma, indicative of volatile depletion. A catas¬ 
trophic mass loss rate of ~10^ kg s“^ around perihelion was seen. Ground-based 
Xingming Observatory spotted the post-perihelion debris cloud. Our morpholog¬ 
ical simulations of post-perihelion images hnd newly released dust grains of size 
a > 10 /im in radius, however, a temporal increase in Omin was also witnessed, 
possibly due to swift dispersions of smaller grains swept away by radiation forces 
without replenishment. Together with the fading prohle of the light curve, a 
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power law dust size distribution with index 7 = 3.2 ±0.1 is derived. We detected 
no active remaining cometary nuclei over ~ 0.1 km in radius in post-perihelion 
images acquired at Lowell Observatory. Applying radial non-gravitational pa¬ 
rameter, Ai = (1.209 ± 0.118) X 10“® AU day“^, from an isothermal water-ice 
sublimation model to the SOHO astrometry signihcantly reduces residuals and 
sinusoidal trends in the orbit determination. The nucleus mass r\j 10 ^- 10 ^ kg, 
and the radius ~50-150 m (bulk density pd = 0.4 g cm“^ assumed) before the 
disintegration are deduced from the photometric data; consistent results were 
determined from the non-gravitational effects. 

Subject headings: comets: general - comets: individual (C/2015 D1 (SOHO)) - 
methods: data analysis 


1. INTRODUCTION 

Near-Sun comets offer valuable chances for studying the physical properties of cometary 
nuclei because they are prone to disintegration due to the proximity to the Sun around 
perihelion (e.g., Sekanina & Chodas 2005). Nearly 3,000 of near-Sun comets have been 
discovered since the operation of the joint ESA/NASA Solar and Heliospheric Observatory 
(SOHO) in 1996, but the majority appeared dim and failed to achieve a high signal-to-noise 
ratio in the data. In many cases observations from SOHO alone contain large astrometric and 
photometric uncertainties. Many questions yet remain unsettled. For example, what are the 
differences between near-Sun comets and other comets in terms of material compositions? 
Where did they originate and how did they evolve into near-Sun orbits? Thus, ground- 
based observations of near-Sun comets are highly encouraged, yet frequently hampered by 
the extremely harsh observational geometry. Consequently it is very rare that a near-Sun 
comet is observed from ground observatories, evidenced by merely £v^ successful examples 
since the start of 21®* century, viz. C/2011 W3 (Lovmoy) (e.g., Lovejoy & Williams 2011, 
Sekanina & Chodas 2012, etc.), C/2012 E2 (SWAN)0, C/2012 SI (ISON) (e.g., Novski et 


*We dismiss the case of C/2008 01 (SOHO), which was serendipitously detected in images of a total 
solar eclipse, after a search based upon SOHO data (Pasachoff et al. 2009). 

^Through private communications, we see that T. Lovejoy managed to obtain 3 astrometric positions 
from his images taken in strong dusk twilight on UT 2012 March 10.38, from Australia. However, his report 
remains largely unnoticed. 
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al. 2012), 322P/1999 R1 and the newly discovered comet C/2015 D1 (SOHO) 

(hereafter 2015 Dl). 

2015 Dl was discovered by W. Boonplod in SOHO's Large Angle Spectrometric Coron- 
agraph (LASCO) images from UT 2015 February 18 (Battams & Knight 2015). Our orbital 
solution (discussed in Section 13.3p to the LASCO astrometry conhrms that 2015 Dl does 
not belong to the Kreutz sungrazing family (Kreutz 1888), or any of the known sunskirting 
families including the Meyer, Marsden and Kracht groups (c.f. Sekanina & Chodas 2005 and 
citations therein). We searched for small bodies with similar orbits via the JPL Small-Body 
Database Search Engine, yet found nothing. The small perihelion distance, q = 6.06 Rq 
(1 Rq = 0.00465 AU), is substantially greater than perihelia of the Kreutz sungrazing fam¬ 
ily ( 5^2 Rq), and close to, but somewhat smaller than, perihelia of the Meyer, Marsden 
and Kracht sunskirting families (mean perihelia all have q ^ 7.7 Rq). Knight & Walsh 
(2013) discriminates near-Sun comets subjected to tidal fragmentation events as sungrazing 
comets. Although 2015 Dl apparently disrupted, the disruption was unlikely tidally driven 
due to the relatively large heliocentric distance. In this manner, we address 2015 Dl as a 
sunskirting comet, rather than a sungrazing one. 

Given the classihcation, we realize that 2015 Dl is a unique sunskirting comet in that it 
is the brightest and the hrst sunskirting comet which was observed from the ground over the 
past half centurj0. This paper presents our photometric, morphological and orbital analysis 
of 2015 Dl. 


2. Observations 
2.1. SOHO 

The SOHO spacecraft is located around the LI point of the Sun-Earth system. The 
LASCO instrument onboard consists of three coronagraphs. Cl, C2 and C3. Only the C2 
and C3 cameras observed 2015 Dl. The C2 and C3 coronagraphs, externally occulted, have 
annular helds of view (FOV) of 1.5-6.0 Rq and 3.7-30 Rq, respectively (Brueckner et al. 
1995). Each instrument is equipped with a hlter wheel, a polarizer wheel, a shutter, and 
a 1024 X 1024 pixel CCD with a pixel scale of 11".9 pixel"^ for C2, and 56".l pixel"^ for 


^Published in Minor Planet Electronic Circular (MPEG) 2015-K84, 
http://www.minorplanetcenter.net/mpec/K15/K15K84.html, 

^Successful ground observations of the other sunskirting comet 322P/SOHO were made on UT 2015 May 
22, later than observations of 2015 Dl. 
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C3. The synoptic C2 data are taken through an orange hlter with bandpass 5400-6400 A, 
whereas the C3 observations are mainly made with a clear hlter with bandpass ~4000-8500 
A. Other hlters are used much less frequently, generally once per day, and these images have 
half resolution (512 x 512). Each camera carries a polarizer wheel having polarizer positions 
of —60°, 0°, and +60°, and takes polarization sequences 1-2 times per day. 

2015 D1 was observed by C3 from UT 2015 February 18.0-21.8, mostly through the 
clear hlter. The C2 camera also continuously monitored it around perihelion from UT 2015 
February 19.6-19.9, all through the orange hlter. Other available data include a few of C3 
blue and orange hlter images, and four triplets of polarized orange hlter images. All the 
LASCO images were processed in a similar way as described in Knight et ah (2010) by use 
of SolarSoftWare (SSW) and SolarSoftWare DataBase (SSWDB)@ in IDL. 

The observational geometry of 2015 D1 from SOHO is illustrated by Figured] 


2.2. Ground-Based Observations 
2.2.1. Xingming Observatory 

We conducted post-perihelion observations of 2015 D1 on UT 2015 March 4, 8, 9 and 15 
via the 10.6-cm f/5.0 refractor attached with an Apogee U16M 4096 x 4096 CCD through a 
photometric standard V-band hlter as part of the Comet Search Program (CSP) of Xingming 
Observatory. The images have a square FOV 4°.0 x 4°.0, and a pixel scale of 3".53 pixeU^. 
Exposures of data taken from the hrst two nights were 60 s and 120 s in duration, whereas 
data from the last two nights had exposures of 120 s only. The image quality varied from 
night to night, generally ~9"-10" FWHM (Full Width Half Maximum). All the observation 
sessions were started from dusk, because of the small solar elongation. Images were hrst 
fully calibrated by subtracting bias and dark current, then were divided by hat-held frames, 
and hnally were normalized by exposure times. The image sequence from each night was 
registered on held stars, and then was shifted following the motion of 2015 Dl. Normalization 
of the sky background was then performed. Finally, the images were median co-added into 
a single frame. We are able to detect an enormous cigar-shaped nebulosity with its west tip 
within ~5' of the predicted positions either given by JPL HORIZONS or our orbit solutions 
regardless of including non-gravitational parameters. It had a dimension of ~ 1° x 0.2°, 
directed approximately east to west. The cloud appeared the most obvious on March 4, 
even discernible in individual frames, and the dimmest yet still sufhcient for visual detection 


^SSW and SSWDB are both parts of the SolarSoft system, http://www.lmsal.com/solarsoft/. 



in the final stacked image from March 15. 


2.2.2. Lowell Observatory 

We attempted to recover 2015 D1 nsing Lowell Observatory’s 4.3-m Discovery Channel 
Telescope (DCT) on UT 2015 March 5. We nsed the Large Monolithic Imager (LMI), which 
has a FOV 12'.3 on a side and a 6.IK x 6.IK e2v CCD. Images were binned on chip 2x2, 
resulting in a pixel scale of 0".24 pixeD^. We obtained three 30 s images using the broadband 
Cousins R filter. Images were trailed at the comet’s rate of motion, and the pointing was 
determined by Lowell’s ephemeris calculator from the orbital solution published by the Minor 
Planet Center (MPEC 2015-D73). LMI was not scheduled to be used on this night, so these 
were the only three images obtained. We removed the bias and applied a fiat field correction 
using images from 2015 February 25, which was the closest night to our observations on 
which science data were obtained. Observing conditions were poor because this was the first 
night following a series of winter storms, so atmospheric seeing was significantly worse than 
normal. Due to the necessity of acquiring images as early as possible following twilight, 
the default focus values were used, so the instrumental point spread function was likely 
suboptimal. 

On the same night we also imaged the comet’s field with Lowell Observatory’s 31-in 
(0.8-m) telescope. The 31-in has a 2K x 2K e2v CCD42-40 chip with a FOV 15'.7 on a side 
and a pixel scale of 0".46. We obtained ten 30 s images with the Cousins R filter trailed 
at the comet’s rate. The bias was removed and the images were fiat-fielded in the standard 
manner. 

Despite that 2015 D1 had a large 3a position uncertainty of ~4' during the observations, 
both of the FOVs are large enough to encompass the region. We visually searched both sets 
of images using several methods but did not find any evidence of the comet. We could 
detect field stars in DCT images to SDSS r magnitude of ~20.0 (Ahn et ah 2012), and likely 
could have detected the comet to magnitude ~21 despite the poor seeing since it would 
have been stationary while the stars were visibly trailed (~3".5 or ~15 pixels). We could 
detect field stars to an SDSS r magnitude of ~19.0 in the 31-in images. This is likely the 
limiting magnitude for any comet non-detection with the 31-in since the stars did not appear 
significantly trailed due to the considerably worse seeing than on DCT. 

Observation condition details from Xingming and Lowell Observatories are summarized 
in Table [H 
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3. Results 
3.1. Photometry 

Although the ground-based observations provided much better resolution than did the 
LASCO C2/C3 cameras, it is impossible to perform photometric measurements with them 
because of the extreme fuzziness of the debris cloud of 2015 D1 as well as the tremendous 
area it occupied. We only conducted aperture photometry of the comet in C2/C3 images by 
using packages in the IDL Astronomy User’s Library (Landsman 1993). 

Apparent magnitudes were converted from the measured fluxes with zero-points of the 
LASCO images, which were calculated based upon trespassing field stars in LASCO data. 
Due to degradation effects of the LASCO detectors, the zero-points have changed slightly 
year by year. Only degradation influences upon C2 orange and C3 clear data have been 
examined exhaustively, as data for these Liters is the most abundant. We cannot hnd out 
any detailed information about the changes in the zero-points of other filters. We adopted the 
temporal zero-point computed by Gardes et al. (2013) to calculate the C2 orange magnitude, 
and the zero-point by Lamy et al. (2013) to calculate the C3 clear magnitude. For other 
Liters we used values given by Llebaria et al. (2006), Knight (2008) and citations therein, 
and further included an uncertainty of ±0.05 mag in error estimates due to the unavailable 
temporal evolutions in the zero-points. 

Because of the low spatial resolutions of the C2/C3 cameras, we used a Lxed angular 
sized aperture, which allows direct comparison of our results to previous studies of near-Sun 
comets (e.g., Biesecker et al. 2002, Knight et al. 2010, Lamy et al. 2013). A circular 
aperture of radius 5 pixels (1^0) was selected for full resolution 1024 x 1024 pixel C2 images 
and 3 pixels (2'.8) was selected for 1024 x 1024 pixel C3 images, to enclose the signal of 
2015 Dl, but at the same time to minimize contaminants from sky background as much as 
possible. Half resolution 512 x 512 pixel images had half sized apertures. The comet appeared 
overexposed in 14 C2 orange Lltered images from UT 18:00-20:48, 2015 February 19, and 
hence we applied saturation corrections, developed by Knight et al. (2012). This likely still 
underestimates the total brightness slightly, but is much closer to the actual brightness. We 
estimate the uncertainties from the saturation correction at < 0.1 mag and are systematic, 
e.g., nearby points have nearly identical saturation correction uncertainties. 

We converted apparent magnitudes my into heliocentric magnitudes Hy, by normalizing 
the distance between SOHO and 2015 Dl to A = 1 AU, and correcting for the phase eLect: 

Hy (rh) = my (rh. A, a) - 5 log A ± 2.5 log [0 (a)], (1) 

where rh is the heliocentric distance, a is the phase angle, and the phase function 0 [a) is 


given by Marcus (2007) as 



1 + fi'f + 25 'f cos a 
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1 + gl + 2g^ cos a 
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Here dgo is the ratio of the dust-to-gas intensity observed at a = 90°, with (Jgo = 1 for normal 
comets, 0 < /c < 1 is the partitioning coefficient between the forward- and back-scattering, 
and gt > 0 and (?b < 0 are respectively the forward- and back-scattering asymmetry factors. 
Marcus (2007) suggested k = 0.95, gf = 0.9 and gb = —0.6 according to observations of 
six comets. The model has been applied widely in recent works regarding observations of 
near-Sun comets, e.g., Knight et al. (2010), Knight & Battams (2014), etc. Although Li 
& Jewitt (2015) found slightly different parameters for C/2010 XI (Elenin), it has minimal 
effects with corrections always < 0.3 mag to 2015 Dl’s data and does not meaningfully 
alter the light curve shape. Several of the comets examined by Marcus (2007) have perihelia 
considerably smaller than that of C/2010 XI and thus the results are likely more comparable 
to 2015 Dl. Therefore, we follow the suggested parameters by Marcus (2007). 

We assign hgo = 1.0 for C3 clear hlter, 5go = 0.39 for C2/C3 orange Liters, and 5go = 10 
for C3 blue hlter, from analysis of 2015 Dl’s color (see Section 14731 for details). Since the 
comet did not experienced strong forward- or back-scattering effects (a ~50-115°), its phase 
function is relatively hat, so the exact choice of hgg always has corrections < 0.2 mag. The 
general shape of 2015 Dl’s light curve would not have been ahected by the phase function 


profoundly. 


The resulting light curve of 2015 Dl is shown in Figure [2^. The comet steadily 
brightened from the beginning at Hy ~ 9, until UT 2015 February 19.4 (denoted as 
At = t — tp ~ —8 hrs, where t is observation epoch, and tp is the perihelion time of 
the comet, UT 2015 February 19.75), when it apparently faded by ~1 mag in ~7 hrs, fol¬ 
lowed by a drastic surge in its brightness to Hy ~ 1.5 through the clear hlter, in ~5 hrs. 
Post-perihelion witnessed a decline in its brightness. The comet was then obstructed by the 
pylon of the coronagraph for 3.6 hrs. Starting from At ~ -|-0.6 day (UT 2015 February 20.3 
or DOY ~ 51.3) it dimmed smoothly on the way out of C3’s FOV. 

Figure [2 ]d shows Hy as a function of rp. We can see that the post-perihelion brightness 
was consistently brighter than the pre-perihelion brightness at the same heliocentric distance, 
by ^ 1-5 mag. The pre-perihelion brightening at rp > 13 i? 0 , oc , was steeper than 
the post-perihelion fading at the same range, oc A turnover point in the inbound 

leg at Tb ~ 13 i ?0 is noticed, where the brightening slowed down to oc The second 

turnover point in the inbound leg occurred at Uh ~ 8 Rq, after which the comet faded despite 
continuing to approach the Sun. Then the hareup took place around perihelion at 6 i? 0 . 









and subsided at outbound r-^ — 6.7 Rq. Similar light curves have been found amongst some 
of the Kreutz group comets (e.g., Knight et al. 2010). Starting from rh ~ 13 Rq in the 
outbound leg, the comet faded steadily toward the end of the LASCO observation. The two 
respectively inbound and outbound turnover points at ~ 13 Rq are very similar to those 
of Kreutz sungrazing comets, which are believed to be related to sublimation of olivines (e.g., 
Kimura et al. 2002). 

As shown in Figure[21 the color of 2015 D1 was initially distinctly different from the color 
of the Sun, yet eventually evolved towards it, indicating that the coma became increasingly 
dusty. We think that this was due to depletion of sodium, which emitted strongly at the 
beginning of the LASCO observation, and faded out gradually. The comet had mean color 
indices Clear — Orange = +0.6, and Clear — Blue = —0.7 (see Section IT75]) . 

Four triplets of LASCO polarizer observations of 2015 D1 are available, however, they 
contain large uncertainties and the number of data points is too small, so the result will not 
be presented. 


3.2. Morphology 

Using SOHO and ground observations, we identihed hve stages in the evolution of 2015 
Dl; 

1. From discovery to At ~ —8 hr (UT 2015 February 18.0-19.4): the comet, initially 
almost stellar, was trailed by a developing faint tail as it brightened. It was similar 
to comparably bright Kreutz sungrazing comets at similar heliocentric distances (see 
Figure H^. 

2. From —8 hrs ^ At < 0 (UT 2015 February 19.4-19.7): the tail weakened and disap¬ 
peared, whereby the comet became completely stellar (see Figure |Dd). 

3. Within 0 < At < +1 day (UT 2015 February 19.7-20.8): the comet developed a new 
tail, which was much more prominent than the pre-perihelion tail, at the same time 
that it brightened by about ~3 mag (Figure [2^). The optocentric region remained 
tight (see Figure St). 


®Note that the time of Figure Hk taken is not within this stage. However, this is the best image which 
shows the existence of the pre-perihelion tail, although it started to weaken. 
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4. From At ~ +1 day until the departure from SOHO^s FOV (UT 2015 February 20.8- 
21.7): the comet maintained its tail, but the optocentric region appeared elongated 
(Figure mi), reminiscent of some notable comet disintegration events such as C/1999 
S4 (LINEAR) and C/2012 SI (ISON) (e.g.. Weaver et ah 2001, Knight & Battams 

2014) . 

5. A week or more after the perihelion passage: multiple ground-based observers reported 
a nebulous cigar-shape object near the nominal position of 2015 D1 (e.g., Masek et ah 

2015) . From these images, we identihed no clear central condensation or a nucleus. The 
object dissipated rapidly as time went by. To our knowledge, no successful observation 
has been reported after mid-March. 

We interpret the physical evolution of 2015 D1 as follows. During its pre-perihelion 
phase, the comet behaved like a typical comet, with ongoing activity producing a dust tail 
(Figure H^). But close to perihelion, dust began to sublimate faster than it was replenished 
and the tail disappeared (Figure Hb)- The rapid brightening around perihelion and the 
subsequent development of a new tail seemed to indicate a sudden surge in activity of the 
comet. Considering the signs of nucleus disintegration depicted in subsequent images, it is 
apparent that such dramatic change of morphology reflects a catastrophic event experienced 
by the nucleus. The low spatial resolutions of SOHO images hamper us from immediately 
looking into details of the disintegration, but it appears that the time from the flareup to 
the ultimate disruption of the nucleus took no more than 1 day (see Section ITT)) . Generally 
speaking, the morphological evolution of 2015 D1 carries many similarities to that of Kreutz 
sungrazing comet C/2011 W3 (Lovejoy) (Sekanina & Chodas 2012). 


3.3. Orbital Determination and Non-Gravitational Effect 

We only used the SOHO astrometric data for orbit determination. 2015 D1 appeared 
too diffuse in ground-based observations, in spite of much better resolutions. Without a 
central condensation it is impossible to conduct astrometric measurements from these data. 

SOHO astrometric measurements were recorded in custom software operated in IDL. 
The basic procedure was to manually select the optocenter of the comet and then allow 
the software to automatically calculate centroids on the 25 closest stars to the comet. This 
process occurred for every image in which the comet was visible. In the case of LASCO 
C2 there were not always 25 stars available, and thus as many as possible were recorded. 
LASCO C3 always has many more than 25 stars available. The limit of 25 stars has been 
selected as an optimum number based on computations of S'CRO-discovered comets in the 
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early part of the SOHO mission. All object locations were recorded at the sub-pixel level and 
passed to an implemented version of the Charon algorithn:0, which reduced the observations 
to a standard MFC format^. 


We used EXORBo to determine 2015 Dl’s orbit. Perturbations by all the eight planets, 
Pluto, and the three most massive asteroids, Ceres, Vesta and Pallas, are included in the 
computation using DE406 ephemerides, although they have basically no influence on solu¬ 
tions. Different weightings were assigned to the observations according to pixel scales. We 
hltered out 10 data points with residuals > 50" as a cutoff. The remaining 412 observations 
all satisfy the residual threshold regardless of considering non-gravitational effects. 


The arc covered by the SOHO observation was larger than any other ^ORO-discovered 
comets so we treated its eccentricity e as one of the free parameters to be solved. We found 
that including the solving of non-gravitational parameters Aj {j = 1, 2, 3), which are defined 
in Marsden et ah (1973)0 from an isothermal water-ice sublimation model, signihcantly 
reduces the sinusoidal trends in astrometric residuals (the differences between the observed 
and calculated positions, a.k.a. O—C residuals, see Figure^]). The trends are irrelevant to the 
selection of the astrometric data, in that filtering further more different sets of measurements 
to stricter residual thresholds or the otherwise, or removing data points apparently close to 
the edge of unblocked regions where diffraction by the occulter and the pylon of LASCO 
might take place do not alter the trend whatsoever. Other factors including infrequent 
resets of SOHO’s onboard spacecraft clock and potential position errors of the spacecraft 
have been fully ruled out. We thus conclude that the residuals are authentic. 


We also found that solving the radial component Ai alone reduces the root-mean-square 
(RMS) of the best ht most noticeably, from ±13".37 to ±10".44. We obtained 6 A 2 > 1 
and 6A3 = 0.31, the relative errors of A2 and A3 respectively, signihcantly larger than 
5Ai = 0.09. Taking into account the poor spatial resolutions of LASCO images, we solved 
Ai = (±1.209 ± 0.118) X 10“® AU day“^ only and simply assigned Aj = 0 for j 7 ^ 1. 
Similarly, JPL HORIZONS gives Ai = (±1.250 ± 0.097) x lO”® AU day^ with Aj = 0 


±ttp://www.projectpluto.com/charon.htm, 

^http://www.minorplanetcenter.net/iau/inf o/ObsFormat.html 

®EXORB, a companion package of SOLEX, is an orbit determination program written by A. Vitagliano, 
available at http://chemistry.unina.it/~alvitagl/solex/ 

^^Marsden et al. (1973) denoted the non-gravitational parameters by Aj (j = 1, 2,3). We feel it necessary 
to change the symbol a little bit so as to avoid potential ambiguity with the used letter Ap, for albedo. 
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for j 7 ^ 1 assumed ^^1. Different weightings and the number of observations hltered by JPL 
HORIZONS may account for the different values. 

Also tested was the forsterite sublimation model by Sekanina & Kracht (2015). But 
we do not prefer that it was the mechanism responsible for the non-gravitational effect 
experienced by the comet, therefore, we did not apply it for the orbit determination (see 
Section 14.61 for details). 

Our solutions to the orbital elements of 2015 D1 are listed in Table [H 


4. Discussion 

4.1. Search for Potential Pre-discovery Data 

We investigated whether serendipitous imaging of 2015 D1 may have occurred. We ap¬ 
plied EXORB to perform multiple iterations of Monte Carlo runs, based upon the random ex¬ 
clusion of a stochastically varying fraction (between 30 - 70%) of the SOHO astrometric data, 
whereby 352, an arbitrary number, Monte Carlo clones of orbital elements of the comet were 
generated. Techniques documented in Clark (2010) were then applied to search for serendipi¬ 
tous pre-discovery imaging of the comet. Using the online Canadian Astronomy Data Centre 
Telescope Products (Gwyn et al. 2012) and the Minor Planet Center Sky Coverage Pointing 
Data dataset, over 600,000 archival images were considered from Canada-France-Hawaii Tele¬ 
scope (CFHT), Panoramic Survey Telescope and Rapid Response System (Pan-STARRS), 
Lincoln Near-Earth Asteroid Research (LINEAR), Spacewatch, Catalina, and Mount Lem¬ 
mon, and approximately 50 smaller surveys, as well as visible and infrared images from 
spacecraft, including Wide-field Infrared Survey Explorer (WISE), and its asteroid-hunting 
portion, NEOWISE, Hubble Space Telescope's Wide Field Planetary Camera 2 (WFPC2) 
and Wide Field Camera 3 (WFC3) catalogues. 1,000 probability clones were generated and 
distributed over a spatial volume consistent with observational errors. Unfortunately no 
recent serendipitous image of any of these clones was identihed. The MPC Sky Coverage 
Pointing Data dataset list three images from Pan-STARRS dated 2014 April 09 and 10 en¬ 
compassing the position of 2015 Dl, however the comet would have been too dim to be 
detected at these early dates, when rh > 5 AU. 
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4.2. Pre-Perihelion Dip in Light Curve 

The SOHO observations showed a dip in the light curve starting ~8 hrs before perihelion 
(Figure [2^). As sodium emission probably contributes signihcantly to the total brightness at 
small Th, one may question whether the dip was caused by the Swings effect (Swings 1941), 
i.e., temporal variation in intensity of cometary emission lines coincident with Fraunhofer 
lines due to Doppler shift. Assuming the entire gas emission was dominated by sodium, we 
use the following equation 

Hv,gas = Hy + 2.5 log [(1 + ^go) Q (rh)] , (3) 

where Q is the normalized g-factor for sodium due to the Swings effect as a function of 
heliocentric radial speed rh, to examine if the dip would be largely removed. We extracted 
a g-factor from Figure 2 in Watanabe et ah (2003) and normalized it to large rh- The 
normalized g-factor has Q = 1 for bh > 80 km s“^ and Q = 0.05 at perihelion. 

We hnd that not only does Equation ([3]) fail to remove the dip in the light curve, but 
also artihcially creates a sharp brightening spike at perihelion. Despite the scatter in the 
data around the bottom of the dip, the maximum dimming ~1 hr prior to the perihelion can 
still be recognized. It is extremely unlikely that perihelion is off by the ~1 hr that would be 
necessary to reconcile it with the Swings effect since this would be more than two orders of 
magnitude larger than the la uncertainty in perihelion (see Table [2]). Moreover, were it due 
to the Swings effect, the dip should occur much more abruptly such that a much sharper 
valley would be formed. More evidence which can help exclude the possibility of the Swings 
effect is that it would shrink the clear — orange magnitude difference centering about the 
dip, which was not seen whatsoever. Therefore the Swings effect is unlikely to be relevant 
to the formation of the dip. 

Likewise, we do not feel that instrumental vignetting can account for the observed dip. 
Admittedly, the minimum of the dip observed by LASCO C3 took place almost exactly 
when the vignetting is locally highest so any inappropriate vignetting correction may cause 
some effects. However, the C2 vignetting is small and relatively constant during the time 
of the cometary light curve dip. LASCO’s vignetting functions are well-established as part 
of the instrument’s calibrations, and accordingly we see no impact on our measurements or 
results that may arise from this correction. We therefore reject the possibility of the LASCO 
vignetting as the reason for the dip; intrinsic activity of 2015 D1 is more likely to be the 
cause. 

We notice that the onset of the dip occurred at rh ~ 8 Rq, following a mild turnover at 
Th ~ 13 i? 0 , consistent with the light curves of Kreutz sungrazing comets. It is thus possibly 
analogous to the turnover in pre-perihelion brightness of the Kreutz sungrazing comets (e.g.. 
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Biesecker et al. 2002, Knight et al. 2010), which is believed to be correlated with the onset 
of sublimation of olivines (e.g., Kimura et ah 2002). The disappearance of the tail around 
perihelion also lends strong support to this idea. 


4.3. Color 


We investigate the color of 2015 D1 based upon hlter magnitude differences. Clear 
magnitudes were determined at the time of orange/blue images by least squares interpolation 
between the nearest clear measurements. The clear magnitude errors were estimated from 
the neighboring clear magnitude errors and were combined with the orange/blue magnitude 
errors using standard error propagation techniques to give a total magnitude uncertainty 
on the color. Figure [3] shows magnitude differences as functions of time and heliocentric 
distances. The non-zero clear — blue and clear — orange magnitude differences, particularly 
pre-perihelion, suggest that the color of 2015 D1 was distinctly different from the color of 
the Sun. However, the color was generally approaching to the solar color gradually as time 
evolved, despite some scatter around perihelion. 


We hrst examine whether the color can be attributed to thermal emission. We approx¬ 
imate dust grains as greybodies. Hence the effective temperature Teg is given by 

1 
4 

(4) 


Tfiff — Cq 


(1 - ^p) 


dear? 


in which Cs is superheat, Ap = 0.04 is a nominal albedo for cometary nuclei! (e.g., Lamy et 
al. 2004), and Sq = 1361 W m~^ is the solar constant (Kopp & Lean 2011), e is the effective 
emissivity, assumed to be unity, and a = 5.6704 x 10“® W m“^ K~^ is the Stefan-Boltzmann 
constant. The influence from thermal radiation is evaluated by 

Fth _ rlJ‘IBx{%s)dX 

F,, /THp0(a)Fo,AdA’ ^ ^ 

where Fth is the thermal emission flux, Fsc is the flux due to scattering sunlight, Fq^\ is the 
solar irradiance spectrum observed at 1 AU, Bx is the thermal emission flux from Planck’s 
law, T is the effective transmissivity of a given hltered optical system, and A is wavelength. 
We calculate T for C2/C3 orange, C3 clean and C3 blue Liters based upon the information 
provided on the LASCO calibration pagim. The 1985 Wehrli solar spectruni^ is used in 


^^http://lasco-www.nrl.navy.mil/index.php?p=content/level_l/lascocal_index 
^^http://rredc.nrel.gov/solar/spectra/am0/wehrlil985.new.html 
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this estimate. Unfortunately there are no available data which can constrain Us; it is a 
function entangled with Tefj and dust grain size (Gehrz & Ney 1992). For simplicity, we 
assume Us = 1.2, which is approximately the average of superheat values of all types of 
comets listed in Gehrz & Ney (1992). 

We compute Fth/Fgc observed by LASGO blue, clear and orange hlters. The simulation 
results are illustrated in Figure [6l from which we can see that influence from thermal radiation 
emission was very limited for all LASGO hlters at Uh ^ 10 Rq. Around perihelion, while G3 
clear data would be affected by thermal radiation the most signihcantly, other hlters would 
still receive negligible thermal radiation. However, since dust grains experience sublimation 
at small rh, the actual equilibrium temperature would therefore be lower than what Equation 
()T|) gives (c.f. Kimura et ah 2002). Hence, we think that the ehect due to thermal radiation 
emission can be ignored. 


Intuitively, were the color of 2015 D1 due to thermal radiation entirely, it is expected 
that the comet would appear redder at smaller rh, which is not observed. Furthermore, as 
indicated in Figure 0 the comet should always appear the brightest in G3 clear images, 
mediocre in G2/G3 orange, and the faintest in G3 blue data, obviously contradictory to the 
LASGO observation (Figure |3]). We can therefore conclude that thermal radiation emission 
is not responsible for the observed color of the comet. 


We next investigate how sodium emission will influence the color of 2015 Dl, since this 
effect is prominent when a comet nears the Sun. Similar to the method described in Knight 
et ah (2010), we add a synthetic rectangular sodium flux ^Na with varying intensity to the 
solar spectrum centered at Na D-line A = 5985 A with a hxed width AANa = 10 A. Since 
AANa is very small, TNa, the mean effective transmission around the Na D-line within AANa, 
can be utilized for simplihcation, such that the modeling magnitude difference between filter 
i and j now becomes 


Arriij 


TxaT^Na, iAANa + / TjUQ^AdA A / J ‘TjUQ^AdAAl 

-^Na^ixa, ,AANa + / T,Wo,,dA J V/‘2:*i"0,AdAy'J ’ 


( 6 ) 


whereby we obtain the results shown in Figure [71 Although only the flux due to the solar 
continuum and the flux due to the sodium emission are taken into consideration, the modeled 
magnitude differences can be matched by varying the intensity of the sodium emission. For 
instance, an intensity of sodium emission ~100 times stronger than the solar continuum at 
5985 A corresponds to an apparent magnitude ~0.4 mag brighter in the G2/G3 orange hlters 
and ~0.3 mag fainter in the G3 blue relative to the G3 clear hlter, which was exactly the 
color of 2015 Dl around UT 2015 February 19.9 (DOY = 50.9). We thus think that the 
sodium emission was a plausible mechanism to account for the color of the comet observed 
in LASGO cameras. 
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The magnitude difference between the clear and orange filters decreased as the comet 
approached perihelion, indicating depletion of sodium emission and the coma becoming in¬ 
creasingly more dusty, i.e. (Jgo increasing. Since the corrections from 6 gQ are generally compa¬ 
rable to the uncertainties in the magnitude data, it is not meaningful to apply a temporally 
varying (Jgo {t) to correct for the phase function. Thus we take the mean magnitude differ¬ 
ences to derive (5go) for clear, orange and blue filters with Figure [Tb respectively. For C3 
blue filter we have (5go = 88, however, several typical cometary emission lines, e.g., C 2 , CN, 
etc., would be transmittable through the bandpass and likely lower this value considerably. 
Thus a conservative (5go = 10 is used. 


4.4. Ejection of Dust Grains 


To understand the morphology of the comet as well as the properties of the remaining 
debris cloud, we employed a Monte Carlo dust model similar to the one used in Ye & Hui 
(2014) to generate synthetic images of the comet. During initial tests we noted the unique 
challenges for the case of 2015 Dl. Firstly, the low spatial resolution of LASCO C3 images 
prevents us from obtaining information about the surface brightness profile of the cometary 
tail. Particularly, the pre-perihelion tail stretched < 4 pixels (i.e. < 4') in these images, 
too small for model comparison. By the time the comet appeared in LASCO C2, the tail 
had begun dimming already. Secondly, for ground-based observations, the combined effect 
from the nebulous nature of the remnant, the lack of a central condensation as a reference 
point, and the large uncertainty of the comet’s position (~1' or ~100 pixels in images from 
Xingming) make it very difficult to directly assess the goodness of the model. Therefore, 
we only focus on matching the general shape of the tail/remnant starting from around 
perihelion. Nevertheless, thanks to the small heliocentric distances at which the dust grains 
were released, the often-significant divergences between different sets of parameters made it 
relatively easy to identify implausible solutions. 


The dynamics of cometary dust grains are determined by the (3 parameter, the ratio 
between the solar radiation force and the gravitational force exerted by the Sun, and the 
initial ejection velocity. The ratio /3, dust grain radius a and bulk density of dust pd are 
related by 


/9 = —, 

Pda 


(7) 


where C = 5.95 x 10“^ kg m“^ is a proportionality constant. After trials with various 
parameter valuations, we found that the post-perihelion shape of the tail/remnant was pre¬ 
dominantly controlled by the generation of small dust particles. Hence, in the following, we 
use the dust ejection model by Crifo & Rodionov (1997) and the upper limit of dust size 
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Omax ~ 1 cm. Assuming a typical bulk density pd = 0.4 g cm ^ (e.g., Richardson et al. 
2007), we have /Smin ~ 1-5 x 10“"^. 

An immediate question is the duration of dust ejection: did the nucleus split instan¬ 
taneously (such that the dust ejection ceased shortly after the comet’s flareup), or did the 
disintegration process last for some period of time? We thus consider three scenarios: 

1. All dust grains were impulsively ejected at the start of the flareup at At ~ —1 hr 
(impulsive ejection); 

2. The dust grains were ejected from the start of the flareup to the peak brightness, i.e. 
— 1 < At < -|-3 hrs (short semi-impulsive ejection); and 

3. The dust grains were ejected from the start of the flareup to the time when signs 
of nucleus disintegration were seen in SOHO images, i.e. —1 < At < -|-1 day (long 
semi-impulsive ejection). 

The simulated particles, isotropically released, were generated using both sets of the 
orbital elements in Table [2] during initial tests. A modihed MERCURY6 package (Chambers 
1999) was used to integrate all particles to observation epochs using the Bulirsch-Stoer 
integrator (Bulirsch 1972, Stoer 1972). Radiation forces are included in the code. Also 
included are gravitational perturbations from the eight major planets, although these cast 
no visible influence on modeling 2015 Dl. We then calculated the positions of simulated 
particles with respect to SOHO or the Earth at epochs of interest to produce the shapes of 
the dust ensembles. 

During tests, by visual inspection, no distinction between the modeled shapes from 
different sets of orbital solutions was detected. We think that different orbital solutions 
affect little the morphological analysis, and therefore applied the solution without Ai. We 
tested /Smax from 5 x 10“"^ to 0.5 using a logarithmically varying interval (i.e., steps of 10“^ 
for /Smax ~ 10“"^, steps of 10“^ for /Smax ~ 10“^, etc.). We selected eight SOHO images from 
UT 2015 February 19 20:06 to February 21 15:18, each separated by about 6 hrs (except 
the hrst two images, from UT February 19 20:06 and February 20 07:18 respectively, are 
separated by 11 hrs, as the comet was obstructed by the pylon of the coronagraph), and 
Xingming images from March 4, 8 and 15 (observations from March 9 were dismissed due 
to a bright background star) for model matching. For SOHO data, the synthetic images 
are essentially a set of segments due to the low resolution of SOHO. The goodness of the 
model is therefore assessed by comparing the distance traveled by different sizes of dust to 
the observed length of the tail (Figure [S]). For Xingming data (Figure [9]), the shape of the 
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remnant seems reproducible; however, we notice that the modeled debris cloud is constantly 
~3' southeast of the actual observed cloud, a phenomenon we attribute to imperfect orbit 
determination. Note that this was present no matter which orbital solution (our own, JPL, 
MPC) was used. The positions of the simulated particles were therefore translated ~3' in 
the northwest direction to align the synthetic images to the observations (Figure |9]). 

The ejection duration is constrained by Xingming data, indicating a quasi-impulsive 
ejection of the dust within 0.1 day (see Figure [HD and endorsing the idea that the destruction 
occurred immediately after the flareup. This is consistent with the analysis by Sekanina 
(2015). The SOHO images, suffering from low spatial resolutions, failed to allow a clear 
separation of different ejection durations, although the length of the tail provides a reliable 
constraint to the lower size limit of the optical dust. An increasing trend of Omin is clearly 
noticeable (Figure dU]). The freshest dust grains had amm ~ 10 pm; it increased at a rate of 
hmin ~ 10“^ mm day“^, and stabilized at ~0.5 mm. This may be explained by observational 
bias: the smallest dust grains in the debris cloud quickly dispersed and thus dimmed beyond 
the observation threshold, and the debris cloud was expanding due to solar radiation pressure 
without replenishment of dust particles. Conversely, larger dust grains expanded more slowly 
and remained observable for a longer period of time. 


4.5. Size Estimate 


4 . 5 . 1 . Nucleus Size from Photometry 


We can estimate the nucleus size of 2015 D1 from the SOHO photometric data. The 
heliocentric magnitude of the comet due to the dust grains which reflect sunlight can be 
extracted, similar to Equation ([3]), by the formula 


-^ydust — Hv + 2.5 log 



( 8 ) 


Then the effective cross-section of the comet can be calculated as 




0.4f J/y.duBt— 
Ap 


(9) 


Here, mQy = —26.74 is the apparent V band magnitude of the Sun. We still use Ap = 0.04 
for the dust grains. Figure ITT] shows Cg as a function of time. We assume that the optically 
thin coma is made of spherical dust grains whose radii range from amin to Umax, and that 
they obey a power-law size distribution, dN oc a“'^da, with a constant 7. Then the effective 
nucleus radius can be solved by 
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1/3 


( 10 ) 


'mm 


From the morphological analysis in Section 14.41 we have 10 /im ^ a ^ 1 cm aronnd 
perihelion. We can constrain 7 from the nniform decline in Cg starting from At ~ 0.6 
day (DOY ~ 51.3) nntil the end of the LASCO observation by assnming that the decline 
was completely attributed to faster dispersions of smaller dust grains accelerated by solar 
radiation forces. The relationship between Cg and the dust size distribution is 



where € is an unknown constant that does not affect the calculation. We know Omin {t) in the 


same interval of time from Figure [lOl A best £t to Equation ffTTl) by MPFIT (Markwardt 
2009) yields 7 = 3.16. It is relatively insensitive to Omax and Cg. For instance, changing 


from 5 mm to 10^ m varies 7 from 3.10 to 3.30. We are confident that distributions 


with 7 = 3.2 ±0.1 encompass the likely range of parameter uncertainties. In comparison, 


distributions with 3.5 < 7 < 4.1 have been found for a large number of comets (e.g., Sitko 


et al. 2011), but 7 = 3.2 ± 0.1 is not uncommon (e.g., Fulle 2004). 

Around perihelion. Equation (fTOj) yields i?N — 0.11 ± 0.01 km. Taking into account 
different assumptions about the albedo (e.g., Lamy et al. 2004 estimates an uncertainty of 


±0.017) yields i?N — 0 . 1 lto;o 2 km, which likely encompasses the original nucleus size. 


4-5.2. Nucleus Size from Non-Gravitational Effect 


Whipple (1950) shows that the nucleus mass Mn can be inferred from the non-gravitational 
acceleration as a result of momentum conservation. The composite non-gravitational param¬ 


eter, A = ± ^2 + is connected to the non-gravitational acceleration by 21 (rh) 


Ag(rt), where g (rt) is the dimensionless empirical non-gravitational momentum transfer 
law from an isothermal water-ice sublimation model (Marsden et al. 1973), which we ex¬ 
ploited in determining Ai- However, the actual mechanism the nucleus of 2015 D1 suffered 
might well be too complicated to be described by any simple models. Additionally, evidence 
suggests that sublimation of olivines (e.g., forsterite) has taken place around perihelion. But 
the contribution to the non-gravitational effect is believed to be very limited (see Section 
14.6p . Given the high uncertainties in the astrometric data, we still apply the empirical law 
g (rh) from the water-ice sublimation model. 



We thus have 


( 12 ) 
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where Q (rh) is the production rate of the dominant mass loss material, i.e. water-ice, having 
molecular mass /i (H 2 O) = 18 u, and k is the dimensionless collimation efficiency, with k = 0 
for isotropic emission and k = 1 for perfect-collimated ejection, and v is the outflow speed 
of gas as a function r^, which is ill-dehned when is small. Applying relationships such as 
those given by Delsemme (1982) and Biver et ah (1998) to a near-Sun scenario is probably 
inappropriate. Instead, we approximate v as thermal speed 


Vth (?^h) 


3fcBy (l-Ap)^e 

/i J Aearl 


(13) 


where ks = 1.3806 x 10“^^ J is the Boltzmann constant. Assuming Ap = 0.04 and e = 1, 
for water-ice sublimation. Equation flT3ll is simplihed to Vth (?"h) = 0.62 km s“^, where 

Th is expressed in AU. 

By no means can k be constrained from the observations, and we somewhat arbitrarily 
adopt K, = 0.5. There is no constraint on the gas production rate of 2015 D1 either, therefore 
the empirical law for long-period comets by Sosa & Fernandez (2011) is applied. We use the 
magnitude only due to gas emission calculated by Equation ([3]) with Q = 1. Hence Equation 
(IT^ yields a mean nucleus mass (Mn) ~ (5.1 ± 3.3) x 10® kg, much smaller than the masses of 
long-period comets studied by Sosa & Fernandez (2011), by 4 orders of magnitude. Assuming 
Pd = 0.4 g cm“®, we have its effective nucleus radius i?N = -^OMn/ (dvrpa) ~ 67 ± 15 m. 

However, in Section 14.5.11 the nucleus size is estimated to be O.III 0 Q 2 km in radius. 
Given the same pd, this yields — (1.3-5.5) x 10® kg, an order of magnitude larger than 
the mass derived from the non-gravitational effect. We consider the following reasons. 


1. The estimated size from photometry includes all the constituents within the aperture, 
which occupies a spatial sphere of ~1.2 x 10^ km in radius around 2015 Dl’s optocenter 
around perihelion, not the nucleus alone. Larger dust grains released earlier would stay 
in the aperture much longer than the entire passage in LASCO C3’s FOV. Therefore 
photometric data give a decent estimate about the initial size, whereas the size estimate 
from the non-gravitational effect tends to give the size around perihelion. We thus 
expect the size estimated from photometry to be signihcantly larger than the one from 
the non-gravitational effect. 

2. The non-gravitational effect might come from a different mechanism other than from 
the isothermal water-ice sublimation model. As mentioned earlier, Ai in the orbital so¬ 
lution fails to remove an obvious leap in the astrometric residuals in declination around 
perihelion (Figure |5]). Together with the photometric data and morphological analysis. 
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this broadly agrees that something catastrophic happened to the comet around peri¬ 
helion. So one might expect that the non-gravitational acceleration emerged predom¬ 
inantly around perihelion, and a smooth and continuous model might have deviated 
from the fact. 

3. The empirical law of gas production rate by Sosa & Fernandez (2011) has never been 
examined at small rh and therefore it may be inappropriate to apply to 2015 D1 directly. 
Alternatively, 2015 Dl’s actual production rate might deviate from the empirical law, 
even though it might still hold at small rj, for other near-Sun comets. 

Given the substantial uncertainties associated with the behaviors of near-Sun comets, 
we think that both methods give acceptably consistent size estimates. We are conhdent that 
the nucleus mass of 2015 D1 was ~10®-10® kg before disintegration, e.g., much smaller than 
most comets studied by Earth-based observers near rh ~ 1 AU. Note that we only used the 
clear images because they had less potential sodium contamination than the orange images, 
and were acquired far more frequently than the blue images. 


4-5.3. Constraints on Post-Perihelion Remnant 

Using the more restrictive Lowell non-detection from the slightly smaller DCT FOV, 
we can also estimate the upper limit of an inactive outbound nucleus of 2015 D1 as follows. 
Assuming the comet has solar color, the SDSS r magnitude converts to V magnitude by 
nir = my-O.lO (Smith et ah 2002), yielding an upper limit of my < 20.16. We then derive 
its effective cross-section with Equation ([9]) and use i?N = \fC^^ to determine an upper 
limit to the remaining nucleus size as i?N ^ 0.6 km. It is necessary to point out that we here 
apply the lAU H-G photometric system phase function by Bowell et al. (1989) for a bare 
nucleus to a = 0° with the slope parameter G = 0.15. 

A radius of i?N ^ 0.6 km is not particularly restrictive, considering that we have pre¬ 
viously shown from both photometry and non-gravitational forces that the pre-perihelion 
nucleus size was ^0.1 km. Different assumptions about the albedo, the phase correction 
(e.g., Lagerkvist & Magnusson 1990 shows that the slope parameter can be off from G = 
0.15 by ~±0.1), or the limiting magnitude (we estimated the comet could have been ~1 
mag fainter than the faintest stars due to trailing), still result in i?N ^ 0.28 km in the most 
restrictive case. 

We next consider the upper limit on an active nucleus radius during post-perihelion 
observations. First we estimate the upper limit to water production rate Q based on our 
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limiting magnitude following the empirical correlation found by Sosa & Fernandez (2011), 
whereby we have Q < 3.3 x 10^® molecules s“^. We next estimate the surface area necessary 
to produce this production rate at = 0.596 AU using the methodology of Cowan & 
A’Hearn (1979) and translate this into an effective radius assuming the comet is active over 
a surface area corresponding to the effective cross section. This yields i?N — 24 m or 50 m 
for the subsolar or isothermal cases, respectively. Given the signihcant assumptions that go 
into this estimate, it is probable that any remaining active nucleus was less than 100 m in 
radius. 


4.6. Mass Loss 


We can investigate the mass loss of 2015 D1 from either photometry or the non- 
gravitational effect. Here we hrst examine the nucleus mass loss from photometry by trans¬ 
forming Equation (IlOp into 

To obtain Ce, hrst, smoothing with 10 neighboring data points is performed to the derived 
Ce (shown in Figure [H]). During tests we found that if too few neighboring data points 
are used, artifacts will be formed from the scattered data around the downhill portion of 
the pre-perihelion dip. On the other hand, there is no signihcant improvement if more 
neighboring data points are included. Next, we take the diherence between each time step. 
With Equation ffTTj) . we then obtain a mass loss rate at each time, as shown in Figure [121 

While optical depth ehects would delay the apparent time of mass loss, this should 
set a reasonable time boundary. We can see that the most rapid mass loss rate occurred 
around perihelion, with ~ 10^ kg s“^. This is consistent with Section iTTI that the post¬ 
perihelion tail was formed during this period in a quasi-impulsive manner. It is noteworthy 
that negative mass loss rates do not necessarily rehect genuine variation; they can be better 
explained by particles continuously drifting out of the photometric aperture or sublimating 
without adequate resupply from the nucleus. We obtain the total mass loss around perihelion 
to be ~10® kg, the same order of magnitude as the original nucleus mass estimated previously. 

We can also investigate the mass loss according to the non-gravitational effect. Equation 
flT^ can be transformed into an ordinary differential equation 


Mn (t) 


Mn (t) V (t) 

^ 21(f) 


(15) 
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where Mn = fiQ. The variables are separable and integrable. We can then solve the ratio of 
mass loss to the initial nucleus mass, Su, during an observation interval from to to fobs, by 


Sm = I — exp 


'_A ^(rh (t)) ' 

. Jto . 


(16) 


Here we choose to to be the time when the hrst LASCO observation of 2015 D1 was made, and 
t varies from the beginning to the end of LASCO observation. We hnd that if the empirical 
water-ice sublimation model is correct, 93.2% of the nucleus mass would have been lost by 
perihelion, and 99.6% would have been eroded by the time 2015 D1 exited LASCO C3’s 
FOV. While we agree that mass erosion is predominantly important to 2015 Dl, judging 
from the photometric and morphological analysis, we suspect that this overestimates the 
mass loss by perihelion. 


We then investigate the mass loss due to sublimation of forsterite and follow the same 
procedures described by Sekanina & Kracht (2015) to calculate ^i,for (•^j.for = 0 assumed 
for j % 1). The sodium model is skipped as we believe that the sodium amount is very 
small and therefore a signihcant mass loss due to sublimation of sodium is highly unlikely. 
During tests we found that the sodium sublimation model gives results very similar to those 
by water-ice sublimation. 

We obtain ^i^for = (3.990 ± 0.362) x 10“^^ AU day“^. Comparisons between different 
models are shown in Figure [121 from which we can see that sublimation of forsterite would 
lead to the comet experiencing rocketing mass erosion once it reached a very small heliocentric 
distance of rh 8 Rq. By perihelion, an overwhelmingly large section of the initial mass, 
89.0%, would be eroded, and the comet would devastatingly lose 99.3% of the mass by the 
end of the LASCO observation. 


However, the forsterite model exaggerates the sinusoidal envelope of the residuals in 
declination pre-perihelion (Figure |5t) and slightly worsens the residuals, RMS = ±11".60, in 
comparison to the water-ice model. Most importantly, we hnd that sublimation of forsterite 
alone fails to support the enormous mass loss experienced by 2015 Dl around perihelion. 
To verify this, we apply equations and parameters in Kimura et al. (2002) to estimate Mn 
due to sublimation of forsterite around perihelion. We obtain the unit area mass loss rate as 
5.7 X 10“^ kg s“^ m“^, which is then multiplied by the surface area of the nucleus, yielding 
Mn ~ 7 X 10“^ kg s“^. Although a porous nucleus would increase the surface area, resulting 
in a larger Mn, yet it is still far too small compared to the mass loss of 2015 Dl around 
perihelion. In comparison, given an isothermal nucleus, the mass production rate of water-ice 
around perihelion is ~0.2 kg s“^ m~^. In order to support the observed peak mass loss rate, 
this would require a surface area of ~1 km^, equivalent to a ~0.3 km radius sphere. This 
is order-of-magnitude consistent with our estimate of the nucleus size before disintegration. 
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Thus, the lack of other better models makes the isothermal water-ice sublimation model the 
best choice for approximation. 


4.7. Mechanism of the Disintegration 


We briefly investigate the breakup of 2015 D1 since it is the only sunskirting comet for 
which there is strong observational evidence for fragmentation while it was being observed. 
We hrst consider the possibility that breakup was caused by tidal disruption due to its 
proximity to the Sun. For a non-spinning fluid body, the Roche radius of the Sun is 3.7 Rq 
for a bulk density of pd = 0.4 g cm“^, whereas 2015 D1 started to fragment at rt > 6.06 Rq. 
A comet experiencing tidal disruption at a distance ~1.6 times larger than the Roche radius 
seems farfetched, although we cannot fully rule out the possibility because of unknown factors 
such as the nucleus density and the way the nucleus spins can affect the actual Roche radius 
(Asphaug & Benz 1996, Richardson et al. 1998). 

We next consider effects from thermal fracture. The timescale of heat conduction from 
the surface to the interior of a spherical rocky body is roughly th ~ R^/nes where Kes ~ 
m^ s“^ is the effective thermal diffusivity typical for rocks. For 2015 Dl, we have th ~ 10^ 
yr. 


The core temperature of 2015 Dl can be estimated from conservation of energy 


Tc 


r(i-7ip)^0 

dt 

_4ecr {tp - to) j 

to 


(17) 


We choose tp — to = 100 yr, and thereby obtain Tc — 90 K for the nucleus core. On the 
contrary around perihelion, with much of its nucleus surface devoid of volatiles assumed 
and sublimation of forsterite taken into account, the equilibrium surface temperature was 
~1640 K; a huge temperature gradient of AT ~ 1550 K from the nucleus surface to the 
interior would be formed. To estimate the established thermal stress we set a nominal 
thermal expansion coefficient, ay ~ 10“^-10“® K“^, typical for common rocks, and a Young’s 
modulus Y ~ 10^-10^^ Pa (e.g., Jewitt & Li 2010, Sekanina & Chodas 2012, and citations 
therein). An overwhelmingly huge thermal stress, dth = ayYAT ~ 10®-10® Pa would 
be generated inside its interior, which is an order of magnitude or more larger than typical 
tensile strengths of cometary nuclei (c.f. Prialnik et al. 2004 and citations therein). Thermal 
fracture and cracking were very likely to occur, whereby preexisting subsurface volatiles were 
exposed, disastrously intensifying the outgassing activity. 


Recently, Steckloff et al. (2015) argued that differential stress within the nucleus interior 
due to dynamic sublimation pressure may have been responsible for the breakup of sungrazing 
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comet C/2012 SI (ISON). This mechanism might be plausible for C/2012 SI (ISON) due 
to the fact that its nucleus had withstood strong outgassing activity for a long period of 
time (> 1 yr) before disintegration, however, there is no evidence that 2015 D1 was similarly 
active. Thus, we favor thermal fracture for 2015 Dl, since sublimation stress is likely orders 
of magnitude smaller than the thermal stress built up within the interior. 


Yet it is still unclear whether the explosion of outgassing directly crumbled the nucleus. 
Even if not, the fate of the nucleus was destined not to survive. Torques exerted by large 
mass loss from the nucleus can lead to rotational instability. Observations suggest that 
rotational breakup is a very common fate for comets in the solar system (e.g., Jewitt et ah 
1997). Concentrating around the perihelion passage, we justify this hypothesis by 


A7?n ~ 


Ib^T^’rot^th’ 


(18) 


which is derived in Li & Jewitt (2015). Here kt ~ 10“^-10“^ is a dimensionless coefficient 
of the torque (e.g., Belton et al. 2011, Drahus et ah 2011), and Prot is the rotation period 
of the nucleus. We assume a rotation period Prot ~ 10® s, which is typical for cometary 
nuclei (Samarasinha et al. 2004). The thermal speed around the perihelion is Uth ~ 1 km 
s“^. By substituting other numbers we obtain APn ~ 0.01-1 m. Combined with Equation 
m, we hnd that around perihelion, it would take the nucleus an extremely short time. 
At = APn/Pn ~ 1-100 s, to achieve such a change in the nucleus radius, which means 
that within such a short period of time, mass shedding due to outgassing would change the 
angular momentum by a signihcant factor. We hence see rotational instability as a plausible 
mechanism for a hnal blow to the nucleus by disintegrating it, provided that it survived 
the outgassing explosion. This agrees with Samarasinha & Mueller (2013) that rotational 
disruption is likely the most common cause for splitting of sub-kilometer sized near-Sun 
comets. Note that rotational disruption is not signihcant at large heliocentric distances 
because the outgassing activity is limited. 


5. SUMMARY 

We present analysis of sunskirting comet C/2015 Dl (SOHO) using observations from 
SOHO, Xingming, and Lowell. We conclude: 

1. Non-gravitational effects experienced by this comet were obvious. Solving Ai in the 
orbital solution improves O—C residuals signihcantly and helps remove the sinusoidal 
trends. We hnd Ai = (1.209 ± 0.118) x 10“® AU day“^, based upon the isothermal 
water-ice sublimation model. The non-gravitational acceleration was unlikely due to 
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forsterite sublimation as there is insufficient sublimation to drive the observed mass 
loss around perihelion. 

2. Photometric data and non-gravitational effects consistently suggest the pre-disintegration 
nucleus mass as Mn ~ 10®-10® kg, and the nucleus size as i?N ~ 50-150 m in radius, 
with pd = 0.4 g cm“^ assumed. 

3. The mass loss was predominantly concentrated around its perihelion passage, with the 
most rapid loss as large as ~ 10^ kg s“^. A significant portion of the mass was 
shed during this time interval, comparable to the original nucleus mass. 

4. Morphological simulation of the comet’s post-perihelion tail indicates that it was 
formed around At ~ —1 hr (UT 2015 February 19.7, or DOY ~50.7) within 0.1 
day, in a quasi-impulsive manner, when the comet suffered from the most rapid mass 
loss. The remnant of the debris cloud was morphologically dominated by smaller dust 
particles. The freshest dust grain sizes were a > 10 pm, with pd = 0.4 g cm“^ as¬ 
sumed. An increasing trend in amin was noticed, which is likely due to the smaller dust 
grains being dispersed more quickly without further replenishment of dust and hence 
dimming gradually beyond the detection threshold. We thus derive a power law index 
7 = 3.2 ± 0.1 for the dust size distribution. 

5. We suggest that the fiareup in brightness was likely triggered by excess thermal stress 
built up within the nucleus interior causing an explosive release of material and ex¬ 
posing subsurface volatiles. The outgassing explosion may have crumbled the nucleus. 
Even if not, subsequent rotational instability of the nucleus could easily lead to its 
disintegration. It would only take the nucleus a very short period of time. At ~1-100 
s, to change its angular momentum by a large factor. 

6 . The huge dip in the light curve starting from ~8 hrs prior to perihelion is not due 
to the Swings effect. Mild turnover points at ~ 13 Rq and the more obvious 
one at rh ~ 8 Rq suggest that sublimation of olivines is likely responsible, which is 
directly supported by the disappearance of the pre-perihelion tail around the same 
time. The subsequent rapid brightening resulted from disintegration of its nucleus, 
which drastically increased the effective cross-section area. 

7. The comet had a color distinctly different from the color of the Sun, in particular pre¬ 
perihelion, but gradually evolved to the solar color. Sodium content and not thermal 
emission was the most likely the cause of the color. Depletion of sodium emission led 
to a final color similar to that of the Sun, which implies that the nucleus exhausted its 
volatiles and the coma turned dusty. 
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8. Ground-based observations from Xingming and Lowell revealed no detectable central 
condensation in the debris cloud 13-24 days after perihelion. The post-perihelion non¬ 
detection from Lowell Observatory restricts any remaining active nucleus size to be 
< 0.1 km. 
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Table 1. Viewing Geometry of Ground Observations 


Date (UT) 

TeP 

Th’’ 

A" 

a'’ 


0’ 

PsAng® 

PsAMV’’ 

X’ 

2015-Mar-04 13:13-13:55 

GSP 

0.579 

0.869 

83.8 

35.5 

154.1 

38.8 

206.3 

2.73-4.07 

2015-Mar-05 02:20-02:26 

L31 

0.596 

0.875 

82.4 

36.6 

154.5 

38.9 

206.5 

2.66-2.81 

2015-Mar-05 02:48-02:51 

DGT 

0.596 

0.876 

82.3 

36.6 

154.5 

38.9 

206.5 

3.59-3.67 

2015-Mar-08 13:15-14:26 

GSP 

0.700 

0.928 

73.6 

42.6 

156.4 

39.9 

208.2 

2.08-3.44 

2015-Mar-09 13:37-14:411 

GSP 

0.729 

0.946 

71.3 

44.1 

156.9 

40.3 

208.8 

2.24-3.59 

2015-Mar-15 13:30-15:29 

GSP 

0.893 

1.075 

59.8 

51.0 

159.0 

43.7 

212.6 

1.72-3.37 


^Telescope: GSP = Xingming Observatory’s 0.11-m refractor; DGT = Lowell Observatory’s 4.3-m 
Discovery Ghannel Telescope; L31 = Lowell Observatory’s 31-in (0.8-m) reflector 

’’Heliocentric distance, in AU 

'’Gometocentric distance to the observatory, in AU 

'’Phase angle, in degrees 

®Solar elongation, in degrees 

’True anomaly, in degrees 

^Position angle of the extended Sun-to-comet radius vector in the plane of sky, in degrees 
’’Position angle of the projected negative heliocentric velocity vector, in degrees 
’Air mass, dimensionless 

ISix images taken later than 14:24 UT were partially obstructed, hence discarded. 


Note. — This table is compiled from JPL HORIZONS. We are aware that discrepancies between 
predicted positions increased over time. The worst case is for the GSP observation on 2015 March 15, 
where the JPL ephemeris differs from the one based upon EXORB by ~0°.8. 
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Table 2. Orbital Elements (Reference: Heliocentric Ecliptic J2000.0) 


Element 

Value without A 

la Uncertainty 

Value with A 

la Uncertainty 

Units 

tp" 

2015 Feb 19.74859 

5.3 X 10-5 

2015 Feb 19.74642 

2.2 X 10-^ 

TT 


0.0284511 

3.1 X 10-6 

0.028219 

2.3 X 10-6 

AU 


1.00099 

2.1 X 10-^ 

1.00142 

2.1 X 10-^ 


A 

69.355 

1.8 X 10-2 

69.582 

2.9 X 10-2 

deg 

0" 

95.924 

1.6 X 10-2 

95.897 

1.7 X 10-2 

deg 

a;' 

235.194 

9.4 X 10-3 

235.635 

4.4 X 10-2 

deg 


- 

- 

±1.209 X 10-6 

1.18 X 10-^ 

AU day-2 


^ Time of perihelion passage in Terrestrial Time (TT) 

Perihelion distance 
^ Eccentricity 
Inclination 

® Longitude of ascending node 
^ Argument of perihelion 

® Water-ice sublimation model by Marsden et ah (1973). Only Ai solved, presumably A 2 = As = 0. 

Note. — The RMSs of the orbital solutions without and with Ai are ±13".37 and ±10".44, respec¬ 
tively. Both solutions have epochs on TT 2015 February 18.00483. 






Heliocentric Distance r^, [Rq] 
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Fig. 1.— Observational geometry of 2015 D1 from SOHO's perspective during its transit 
in SOHO's FOV. The vertical dotted line in each panel marks the perihelion time t-p (UT 
2015 February 19.75). 
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(a) 



(b) 

Fig. 2.— V band heliocentric magnitude of 2015 D1 observed by SOHO/LASCO as (a) a 
function of time and (b) a function of heliocentric distance. Point symbols correspond to 
telescopes and points are color coded according to hlters, as shown in the legend. The upper 
panel labels perihelion by a vertical dotted line. The two arrows in the lower panel sketch 
the direction of the comet’s evolution. 
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Fig. 3.— Color of 2015 D1 observed by SOHO/LASCO as (a) a function of time and (b) a 
function of heliocentric distance. The upper panel shows the perihelion moment by a vertical 
dotted line. 
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Fig. 4.— Morphological evolution of 2015 D1 observed by SOHO. The top two panels, (a) & 
(b), are LASCO C2 images and the bottom two, (c) & (d), are LASCO C3. In each panel. 
North is to the top and East to the left. The blue arrows point to the projected negative 
heliocentric velocity vector, and the white arrows point to the projected anti-solar direction. 
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Fig. 5.— Plots of O—C residuals in right ascension and declination as fnnctions of time in 
different orbit determinations. The left panel (a) shows residuals from the pure gravitational 
solntion, the middle one (b) shows residnals from the non-gravitational solntion based on an 
isothermal water-ice snblimation model, and the right one (c) are residuals from the non- 
gravitational solntion with a forsterite snblimation model. A sinusoidal shape in the left 
panel is clearly seen. Althongh signihcant residnals still exist, the solntion with a water-ice 
snblimation model overall gives the best RMS and removes the pecnliar trends presented in 
the left panel. Each panel marks the perihelion of 2015 D1 by a vertical dotted line. Note 
that the three panels have different ordinate scales. 
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Heliocentric Distance [Rq] 


Fig. 6.— Assessment of influences from thermal radiation approached by examining the 
ratios of thermal emission flux to solar continuum flux Fth/Ftc as a function of heliocentric 
distance rh, observed in different SOHO/LASCO bandpasses. The closer to the Sun, the 
more influential thermal radiation is. 
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Fig. 7.— Modeled influences from sodium emission observed in different SOHO/LASCO 
bandpasses. Details are discussed in Section 14.31 Note that C2 orange and C3 orange show 
no obvious differences and therefore overlap each other. 
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Fig. 8.— Termini of dust grains at different /^max in four selected SOHO images: (a) UT 
2015 February 20 07:18; (b) February 20 18:06; (c) February 21 5:42; and (d) February 
21 15:18. The dust models shown here were generated using the impulsive ejection model. 
The difference between impulsive ejection and short/long semi-impulsive ejection is not dis¬ 
tinguishable in SOHO images. Ticks are plotted in the interval of 10', and f3 values are 
indicated on the plots. The images are oriented such that north is up and east is left. 
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(b) 




Fig. 9.— Xingming image on 2015 March 4 overlaid with the best dust models (contours) of 
(a) impulsive ejection (t —tp ~ —1 hr); (b) short semi-impulsive ejection (—1 ^ t —tp < +3 
hrs); and (c) long semi-impulsive ejection (—1 ^ t — tp < +1 day). The results for March 8 
and 15 are largely identical. Dust models are translated ~3' northwest to counter the offset 
presumably introduced by an imperfect ephemeris. The model agrees with the observation for 
the cases of impulsive and short semi-impulsive ejections (i.e. ejection duration < 0.1 day). 
Ticks are plotted in the interval of 10'. The images are oriented so that north is up and east 
is left. 
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Fig. 10.— The temporal decrease of /Smax as seen in SOHO and Xingming data. 
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Day of Year 2015 

Fig. 11.— Temporal variation of effective cross-section area of 2015 D1 against time from the 
LASCO observation. The vertical dotted line labels the perihelion moment. Point symbols 
correspond to telescopes and points are color coded according to filters. 
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Fig. 12.— Mass loss rate calculated from photometric data. Only C3 clear data are used 
because of the adequate number. The perihelion moment is marked by a vertical dotted 
line in the middle of the graphic. Negative values in the plot should not be regarded as the 
authentic mass loss rate of the nucleus, but that the mass loss rate decreased due to particles 
drifting out of the photometric aperture. 
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Fig. 13.— Modeled cumulative mass loss ratios from two different non-gravitational mo¬ 
mentum transfer laws, i.e. water-ice and forsterite sublimation. The models are labeled on 
the plot and detailed discussions are in Section 14.61 











